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Abstract 

Nonlinear effects on the early stage of phase ordering are studied using Adomian's 
decomposition method for the Ginzburg-Landau equation for a nonconserved order 
parameter. While the long-time regime and the linear behavior at short times of 
the theory are well understood, the onset of nonlinearities at short times and the 
breaking of the linear theory at different length scales are less understood. In the 
Adomians decomposition method, the solution is systematically calculated in the 
form of a polynomial expansion for the order parameter, with a time dependence 
given as a series expansion. The method is very accurate for short times, which 
allows to incorporate the short-time dynamics of the nonlinear terms in a analytical 
and controllable way. 
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1 Introduction 



The study of phase ordering kinetics is of fundamental importance to the un- 
derstanding of nonequilibrium dynamics of phase transitions and is of interest 
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in diverse branches of physics, ranging from cosmology and elementary par- 
ticle physics to the different areas of condensed matter physics. The typical 
situation for phase ordering occurs when a system consisting of a two-phase 
mixture in a homogeneous phase is rapidly driven across the critical coexis- 
tence temperature T c into a nonequilibrium state. Fluctuations around the 
initial homogeneous state will develop and the system will ultimately break 
into domains of different phases in space, forming a new equilibrium state. 

The dynamics of the time evolution towards equilibrium is a much-studied 
problem is nonequilibrium statistical mechanics - for reviews see Refs. [1] and 
[2] . If the order parameter that characterizes the phases of the system is not 
conserved, the time evolution of <fi is usually described by the phenomenological 
Ginzburg-Landau (GL) equation [2] 

^ = . M M + CRj) , (!) 



where F[cf>] is the coarse-grained free energy functional, and M is a constant 
mobility, and ((x, t) is a noise term that mimics the thermal fluctuations. For 
the purposes of the present paper, we take for F[(f>] the usual Landau- Ginzburg 
double-well form 



F[</>] = jd*x [i^(V0) 2 + ie0 2 + l 



u i 



(2) 



Here, R 2 is a measure of the interaction range, e = (T — T c )/T c is the reduced 
temperature, and u is a positive constant. In a situation of a temperature 
quench, e = — | e|, the equilibrium (long time) values of the order parameter 
correspond to 2 = | e\/u, the binodal points. The unstable (spinodal) region 
corresponds to 2 < | e|/(3tt). There are several time regimes associated with 
the time evolution towards equilibrium. When <fi at t = 0, one has the 
linear regime. The nonlinear terms in Eq. (1) are small for small t and can be 
neglected; the equation becomes linear [3]. In this case, the equation can be 
integrated very easily and the (noiseless) solution can be written as 

(f)(x, t) = J d 3 k exp [ik ■ x — M(R 2 k 2 + e) f] <j>(k), for t « 0, (3) 



where (f>(k) is the Fourier transform of <p(x,t = 0). In a situation of a tem- 
perature quench, there is an exponential growth of the order parameter for 
wave-lengths corresponding to k 2 < \ e\/R 2 . Because of the fast exponential 
growth for long wave-lengths the nonlinearities rapidly set in, the growth pro- 
cess is slowed down and for longer times interfaces start to be produced, and 
the linear theory cannot be employed anymore. The time scale that character- 
izes the onset of the nonlinear regime depends on M and R 2 . For much longer 
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time scales, one will reach the final stage of the coarsening process and the 
reach of equilibrium. 

In the present paper we are mostly interested in the study of the breakdown 
of the linear approximation at early stage of the phase ordering. In particular, 
we are interested in studying systems with initial configurations that consist 
of domains where the order parameter is not necessarily small. This situa- 
tion occurs for example if the initial state of the system consists of domains 
composed of metastable matter, i.e. \e\/u > 4> 2 > |e|/(3u). Our aim here 
is to provide a semi-analytical method to study the short-time evolution of 
nonlinearities in the Ginzburg-Landau equation for such situations. The mo- 
tivation is of course to obtain an understanding of the time and length scales 
of the breakdown of the linear approximation in a systematic and controllable 
way, without resorting to some sort of discretization methods which can easily 
result in massive numerical computation. 

Specifically, we propose to employ Adomian's decomposition method [4] to 
represent the solution of Eq. (1) in terms of a functional series expansion. The 
method has been used in several instances [5] and more recently it has been 
used for the KdV equation [6] and the coupled system of the Schrodinger-KdV 
equations [7] . For the GL equation the use of the method has been outlined in 
Ref. [8], but no explicit solutions were obtained nor comparisons with complete 
solutions were done. In the present letter we show that the decomposition 
method can be a useful analytical tool to investigate the dynamics of the 
early stage of phase ordering kinetics. We focus primarily on situations were 
the linear theory certainly fails, as in the case when the initial state consists 
of large domains with metastable matter. Since the early works of Ref. [9], the 
study of the early stage of phase ordering has gained renewed interest in recent 
years [10]- [15]. Another direction that the decomposition method can be useful 
is in the study of the scaling behavior at early stages of phase ordering [16], a 
subject of current interest for thin films roughness studies [17]. 

In this first application, to illustrate the applicability and reliability of the 
method for the GL equation, we use a simple double-well free energy func- 
tional and present numerical results for the one-dimensional case only. The 
initial condition is chosen to be an oscillatory function to mimic a configu- 
ration of domains with large order parameter. A more detailed investigation 
at higher dimensions and different free-energies and initial configurations is 
reserved for a future publication. In the next Section we explain the decom- 
position method for the GL equation and obtain the lowest-order terms of 
the series expansion. In Section 3 we show that in the linear approximation 
the series can be explicitly summed and gives the exact result. For the non- 
linear regime we use lowest-order terms of decomposition to compare results 
with the numerical solution of the GL equation. Our conclusions and future 
perspectives are presented in Section 4. 
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2 The decomposition method for the Ginzburg-Landau equation 



For the free energy of the form given in Eq. (2), the Ginzburg-Landau equation 
(we neglect the noise term since in a quench situation e < the thermal 
fluctuations being proportional to y/T become smaller after the quench) can 
be written as 



d(j)(x, t) 
dt 



( 7 V 2 + /3)0(f,t)-«40(f,t)], (4) 



where A[cf>] is the nonlinear term 

A[<f>(x,t)] = <f> 3 (x,t) , (5) 

and 7 = MR 2 ,/3 = M\e\ and a = Mu. Adomian's decomposition method 
consists in expressing in the form of a functional series expansion, where 
the terms of the expansion are determined recursively. Specifically, the itera- 
tive procedure consists in writing (f>(x,t) and the nonlinear term A[<f>] as the 
expansions 

oo oo 

0OM) = S 0»OM) > A [^\ = X]A>[0O,01,---,0n], (6) 

n=0 n=0 

where the A^s are the Adomian polynomials, which can be determined from 
A„[0 O , 01, ■ ■ ■ , <f>n] = ^ | J^[0o + A Vi + • • • + A"0 n ] | . (7) 

Rewriting Eq. (4) as 

t t 

<f>(x, t) = <f>(x, 0) + (7 V 2 + (3) j dt 1 <f>(x, t')-a j dt' A[<f>(x, t')} , (8) 



and substituting the expansions of Eq. (6), one obtains 

00 * 00 

Y, <t>n{x, t) = <f>(x , 0) + ( 7 V 2 + P) / dt' Y <t>n{x, f) 
n=0 q n=0 

i 00 

-ajdt'Y, An[Mx, t'), 4>i(x, *'),•••, 0nOr, t')] ■ (9) 



n=0 
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If one identifies (f>o(x,t) = <p(x,0) = <po{x), Adomian's iteration procedure 
consists in matching the (n+ l)-th order term on the l.h.s. of this equation to 
the n-th order term on the r.h.s. as 



0n+l (X, t ) = ( 7 V 2 + P) J dt' (j) n (x, t') 



t 

-ajdt' A n [(j) (x, t'), 0i (x, *'),•••, <t>n{x, t')} . (10) 



For the case when A[cf>] is a cubic form, the A n 's can be written in a compact 
way as 



1 f d n n 

n 

= E 0(i + 3 + k - n) MAk , (11) 

i,j,k=0 

where 6 is the Heaviside step function. The general expression for the (n+l)-th 
term of the expansion of <j)(x, t) is given as 



<j> n+1 {x , t)= J dt' [(7V 2 + ff)<j> n {x , t') 





a jr ^(^+J + A;-n)0 ^ (f,^')0 J (^,^ / )0A ; (^,^ , )] • (12) 

i,j,k=0 



In general, it is not possible to express the solution in closed form. As will be 
shown in the next section, when the nonlinear term is neglected, the series can 
be summed and the solution is precisely Eq. (3). The solution will be given in 
terms of a power series expansion in t as 

t n 

4> n (x,t) = 4> n {x) — , (13) 



where the n 's are easily obtained from Eq. (12). The first few terms are given 
by 

= (14) 
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1 (f) = ( 7 V 2 + f3) o (x)-a0 o i (f), 



Ix) (pAx) 



(15) 
(16) 
(17) 



(7 V 2 + f3) fa(x)-3a 
3 (f ) = (7 V + P) fa(x) -3a [0 o (f) 2 (f ) + 2 (f ) 2 (f) 
4 (f) = (7 V 2 + (3) fa(x) 

- a [4>l(x) + 30g(f)0 3 (f) + 60 O (f)0l(f)0 2 (f)] , 
5 (f) = (7 V 2 + /3) 04 (X) 

- a 30 2 (f)0 2 (f) + 30o(f)0 4 (f) + 60 O (f)0l(f)03(^) 

6 (f) = ( 7 V 2 + /3) 5 (x) 

30 2 (f)0 3 (f) + 30 2 0! + 30 2 (f)0 5 (f) + 60o(f)0 2 (f)03(f)l • (20) 



(18) 
(19) 



— a 



In the next section we present explicit solutions for the one- dimensional case. 
We obtain the exact result in the linear approximation and for the full non- 
linear equation we compare results obtained with the decomposition method 
with the full numerical solution of the Ginzburg-Landau. 



3 Explicit solutions 



We start considering the linear equation (a = 0). The n-th term is given by 
1 



n (f,t) = - (iV 2 + (3)t o (f). 
n\ 1 J 



(21) 



The series can be explicitly summed, giving 

00 

0(f , t) = ]T n (f , t) = exp [(7V 2 + 0) t] 0o (f). 



n=0 



(22) 



Using the Fourier transform of the initial condition and applying the expo- 
nential operator to it, one obtains precisely Eq. (3) with e = — | e|. 

Next, we consider the nonlinear terms. We use an initial configuration that 
mimics a system with domains of phases characterized by an order parameter 
that is not small. We also consider a quench to zero temperature (e = —1). For 
u 2 = 1, we have that the unstable (spinodal) region corresponds to — yT/3 < 

< +^1/3, and the equilibrium values for the order parameter are (p eq uu = 
±1. Specifically, we use as an initial configuration a periodic function of the 
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form 

0(x, 0) = 0o cos(ax) 



(23) 



In the regions of space in the neighborhood of x — ±nir/a, n = 0, 1, • • • the 
order parameter is close to 0o; as one departs from these regions the order 
parameter starts diminishing until it reaches zero, independently of the value 
of O . Note that the average value of <f>(x, t) is zero over a large region of space. 

It is instructive to take a look at the first few terms of the 0's of Eqs. (14)- (20). 
We present the explicit forms of the n up to n = 4, 



01 (X) 



= 0o cos(ax), 
= 0o cos(ax) 



— a 



sin (ax 

x 



04 (^) 



(24) 
(25) 



(26) 



, i - a 0o cos 2 (ax 
= 0ocos(a:r) {P — a 2 1) — a P — 5 a 2 7) 0q cos 2 (ax) 
+ 2 a 2 0q cos 4 (aa;) — 6 a a 2 7 0q sin 2 (aa;) , 

= o cos(ax) {p — a 2 7) + a 2 (vi (3 — 31 a 2 7^) 0q cos 4 (ax) 

- 9 a 3 0° cos 6 (ax) + 24 a a 2 7 (-,3 + 3 a 2 7) 0q sin 2 (( 
+ a (-9 P 2 + 26 a 2 /3 7 - 33 a 4 7 2 + 58 a 2 a 7 2 sin 2 (a 
x 0q cos 2 (ax) , 
= 0o cos(a x) I (/9 — a 2 7) — 3 a 3 (31 /3 — 67 a 2 7^ 0q cos(a x) 6 
+ 40 a 4 O 8 cos(a x) 8 + a 2 O 4 cos(a x) 4 [71 p 2 - 260 a 2 /3 7 

+ 409 a 4 7 2 - 598 a 2 a 7 0q sin 2 (a x)] + 348 a 2 a 4 7 2 0q sin 4 (a x) 
-6aa 2 7 (l3/? 2 - 66a 2 /?7 + 117 a 4 7 2 ) 2 sin 2 (a:r) 

+ a<pl cos 2 (ax) [-19/? 3 + 83a 2 /3 2 7 - 189a 4 /?7 2 + 253a 6 7 3 

- 2 a a 2 7 (-253 P + 929 a 2 7) 2 sin 2 (ax)]). (28) 



(27) 



The effects of the nonlinearity are encoded by the terms proportional to a and 
involve powers of 0o higher than two. It seems clear that when O is small, 
the higher powers of 0o will give small contributions. We verified explicitly 
that at small times one is still in the linear regime at short times when 0o is 
small and the decomposition method gives the correct answer for all values 
of x. When 0o is large, the nonlinearities will modify the linear term - in the 
regions close to the coordinate points x = ±im/a, n — 0, 1, • • • - at relatively 
low order n , even at very short times, as we shall see shortly. Therefore, 
despite the average value of is zero over a large region of space, for local 
regions with large values of the local approximation breaks down. This is 
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agreement with the discussions in Refs. [11] and [12] on the breakdown of the 
linear approximation at different lenght scales. 

For longer times, it seems that one has to go to large orders in the expansion to 
obtain a good approximation to the solution. The situation can be improved if 
one is able to perform a (partial) summation of the series. In the literature [18] 
Pade approximation has been used with success to improve the convergence 
of the series. We are not going to follow this path here, since we are interested 
in the short-time behavior of the solutions. 

In order to obtain a quantitative estimate of the effects of the nonlinear terms, 
we compare results with a numerical solution of the GL equation. The numer- 
ical solution is obtained using a semi-implicit finite-difference scheme for the 
time evolution and a Fast Fourier Transform for the spatial dependence [19]. 
Specifically, we use a domain of length L, with node points Xj, j — 0, • • • , N, 
with spacing given by h = L/N. The value of <f)(x,t) at point Xj at time t n is 
denoted by 0™. The spatial and time derivatives are discretized as 

d 2 <PM ~ 20," + g-i d<j>(x,t) 0? - 0- 1 m) 

dx 2 h 2 ' dt At ' 1 ' 



Next, we write the Fourier series for 0™ and t/(0") = /30™ — a(0™) 3 as 

$1 = E ex P V^ k ) U W) = E ex P V^ k ) un ^- ( 3 °) 



Substituting these into the GL equation, one can write the following semi- 
implicit equation for at momentum k and time n 



i 

7A fc At 



ln-1 



(k) + At U n ~\k) 



(31) 



where = [2 — 2 cos(2nk /N)]/h 2 . This scheme has been shown to be very 
efficient previously [19], and seems to be adequate for our purposes here. 

In Figure 1 we plot the solutions at short times (t < 1/M) obtained with 
the decomposition method, for the linear (dashed lines) and nonlinear (long- 
dashed lines) approximations, and with the numerical method just described 
(solid lines). We present the results for x — 0. We use a = 7r/4, L = 8 and 
0o = 0.5 and O = 0.75 for the initial condition in Eq. (23); the first value 
corresponds to unstable matter and the the second to metastable matter in 
the region close to x = and x = L. Note that the average value of over 
our spatial length L is zero (we use periodic boundary conditions). 
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^ = 0.5 ^ = 0.75 




0.25 0.5 0.75 1 0.25 0.5 0.75 



Mt Mt 

Fig. 1. The solutions at short times at x = of the noiseless GL equation. The solid 
line represents the numerical solution using Eq. (31) and the dashed and long-dashed 
lines represent the solutions obtained with the decomposition method up to 6-th 
order, in the linear and nonlinear approximations respectively. 

From the plots in the Figure, one clearly sees that the linear approximation 
fails dramatically at all time scales. It is also clear that the decomposition 
method describes very well the exact solution at short times, up to the char- 
acteristic "diffusion" time scale given by 1/M. 

In order to analyze the breakdown of the linear approximation and the rate 
of convergence of the expansion we consider the partial sum 

N 

<t> N {x,t) = $> n (x,t). (32) 

n=0 



The exact solution is obtained for N — > oo. In Table 1 we present the ratios 
of 4>n to the exact solution (f> exa ct(x = 0,t) for the linear (4> l N ) and nonlinear 
(0jv) approximations for t = 1/(2M) and t = 1/M. One clearly sees again 
that the decomposition method provides an excellent approximation to the 
solution at short times at relatively low order. For iV = 4, in the nonlinear 
case, it approximates the solution to 2.5% accuracy at t — 1/(2M) and to 
10% at t — 1/M. The Table also indicates that convergence becomes slower 
after such an accuracy has been achieved and a resummation seems necessary 
to further improve convergence. 
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Table 1 

Ratios of the linear (f) l N (x = 0,i) and nonlinear 4 7 n( x = 0,t) approximations to the 
exact solution 4'exact 

(x = 0,t), for t = 1/(2M) and t = 1/M. 



t = 1/(2M) 



t = 1/M 



N 


tNI 4 s exact 


4% / 4'exact 


YN 1 4 exact 


4^ 1 4 exact 


1 


1.088 


0.997 


1.185 


1.036 


2 


1.178 


1.031 


1.478 


1.146 


3 


1.193 


1.026 


1.576 


1.115 


4 


1.194 


1.025 


1.600 


1.100 


5 


1.195 


1.025 


1.605 


1.100 


6 


1.195 


1.025 


1.606 


1.100 



4 Conclusions and Future Perspectives 



We have proposed to use Adomian's decomposition method to study in an 
analytical way the early stage of phase ordering in the context of the Ginzburg- 
Landau equation for a nonconserved order parameter. The method seems to 
be particularly useful for studying the time evolution of the order parameter 
for systems with regions in space containing metastable matter, where the 
linear approximation breaks down. Although we have used a particular free 
energy functional and a simple initial configuration, they were sufficiently 
representative for the purposes of the present paper. 

The results presented here apply to a model A type of equation, with a double- 
well free-energy functional. We expect the method to be useful for other types 
of phase ordering models, such as model B, governed by a Cahn-Hilliard type 
of equation. Work in this direction is in progress. The coupled problem of a 
conserved and a nonconserved order parameter (model C) seems to be par- 
ticularly interesting, since very little is known about the short time dynamics 
for nontrivial couplings. 
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